arXiv:1504.05049vl [cond-mat.str-el] 20 Apr 2015 


Inexhaustible physics of the helical magnet MnSi: field evolution of the magnetic 
phase transition inferred from ultrasound studies 

A.E. Petrova and S.M. Stisho\0 
Institute for High Pressure Physics of RAS, Troitsk, Russia 

The longitudinal and transverse ultrasound speeds and attenuation were measured in a MnSi 
single crystal in the temperature range of 2 - 40 K and magnetic fields to 7 Tesla. The magnetic 
phase diagram of MnSi in applied magnetic field appears to depend on the experimental setups, 
which is related to a difference in demagnetization factors arising due to the disc shape of the 
sample. 

The magnetic phase transition in MnSi in zero magnetic field is signified by a quasi discontinuity in 
the cii elastic constant, which varies significantly with magnetic field. It is notable that the region 
where the Cn discontinuity almost vanishes closely corresponds to the extent of skyrmion phase 
along the magnetic to paramagnetic transition. This implies that the cn elastic constant is almost 
continuous through the transition from the skyrmion to paramagnetic phases. A recovery of the 
discontinuity of cn and enhanced sound absorption occur at the crossing of the phase transition line 
and the line of minima in cn. The powerful fluctuations at the minima of cn make the mentioned 
crossing point similar to a critical end point, where a second order phase transition meets a first 
order one. 

The skyrmion domain in the case of a ’’perpendicular” setup with a smaller demagnetization 
factor has a reduced temperature range, which suggests that the magnetic field inhomogeneity plays 
an important role in the skyrmion occurrence and, hence, opens a way of skyrmion manipulation. 

The small anisotropy of the shear moduli in the (001) plane found in the parallel setup is most 
probably also caused by the magnetic field inhomogeneity, which distorts the hexagonal symmetry 
of the skyrmion crystal. 


I. INTRODUCTION 

As is well known, the intermetallic compound MnSi 
experiences a phase transition to a spin ordered phase at 
~29 K— The magnetic order in MnSi at zero magnetic 
field is identified as a long period ferromagnetic spiral 
or a helical spin structure with a helical twist caused 
by the Dzyaloshinski-Moria interactiorii^. The helical- 
paramagnetic phase transition at zero magnetic field, 
thought to be of second order, is now recognized as a 
first order transition. Some experimental data in favor 
of a first order transition in MnSi were obtained long ago 
and were subsequently supported by theoretical analy¬ 
sis^. But these indications were mostly ignored, possibly 
because of potential contradiction to the attractive con¬ 
cept of the existence of a tricritical point in MnSi at high 
pressure^. However, new data on the heat capacity, ther¬ 
mal expansion, electrical resistivity, the elastic properties 
and from neutron scattering clearly evidenced the first 
order nature of the phase transition in MnSi^“— . Fig. [T] 
illustrates the behavior of various physical quantities at 
the phase transition in MnSi. Sharp peaks on the low 
temperature side of shallow maxima of fluctuation origin 
(Fig. [Hi) are characteristic for the heat capacity, ther¬ 
mal expansion coefficient, temperature coefficient of the 
electrical resistivity, and the ultrasound attenuation^^— . 
Elastic moduli near the phase transition in MnSi demon¬ 
strate small quasi discontinuities on the low temperature 
side of the fluctuation minima (Fig. [T]d}^. This discon¬ 
tinuity clearly evidences the first order character of the 
magnetic phase transition in MnSi. Of course, a discon¬ 
tinuous change of elastic constants may be observed also 


at some second order phase transitions but that kind of 
transition can not be accompanied by a diverging heat 
capacity as it does in case of MnSi^ii. 

Currently the Brazovski theory of fluctuation-induced 
first order phase transitions has become very popular in 
the analysis of the MnSi syste m^^d^ . 

On application of magnetic fields the helical spin struc¬ 
ture of MnSi transforms to a conical one and finally to a 
spin polarized paramagnetic structure^^.A rough scheme 
of the H—T phase diagram of MnSi is shown in Fig.jH- A 
small pocket of a so-called A-phase exists in the vicinity 
of the transition to the paramagnetic phase at magnetic 
fields ~0.1-0.2 T— . Recent neutron scattering studies 
have revealed features of the skyrmion spin structure of 
the A-phasei^. Heat capacity studies of the magnetic 
phase transition in magnetic fields permit the authors of 
Ref to claim the existence of a tricritical point at ~0.4 
T and ~28 K on the T — H equilibrium line (point 3 in 
Fig. (H), corresponding to the coexistence of the mag¬ 
netic and paramagnetic phase of MnSi. The latest ultra¬ 
sound studies have revealed a rather sharp change of the 
elastic behavior of MnSi at the phase line separating the 
skyrmion domain from the conical phased. However, it 
is still not quite clear what the nature of the phase tran¬ 
sitions from the skyrmion crystal to the paramagnetic 
and conical phases is. Or, equally, what are the charac¬ 
teristics of the crossing points 1,2 (Fig. IH) , where the 
helical/conical phase line enters and exits the skyrmion 
domain. 

The point is that depending on the nature of the phase 
boundaries these crossing points could be simple triple 
points, or multicritical points, or critical end points. The 
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whole situation crucially depends on the nature of the 
phase transition from the skyrmion phase or phase A to 
the paramagnetic phase. Indeed it would not be surpris¬ 
ing if the nature of the indicated phase line differs some¬ 
how from those of both the helical and conical transitions 
to the paramagnetic phase. 

To resolve these issues we performed ultrasound stud¬ 
ies of MnSi in magnetic fields 0 - 7 T and temperature 
2 - 40 K with an emphasis on the magnetic and A-phase 
transitions. 

As discussed above, elastic moduli at the phase tran¬ 
sition in MnSi demonstrate small quasi discontinuities 
on the low temperature side of the fluctuation minima 
(Fig.[T})), which can be used for identifying the character 
of the phase transition^. 



FIG. 1: (Color online) Schematic view of the behavior of a) 
the heat capacity, thermal expansion coefficient, temperature 
coefficient of electrical resistivity, ultrasound attenuation and 
b) elastic moduli at the phase transition in MnSi c) schematic 
magnetic phase diagram of MnSi. 



FIG. 2: (Color online) Temperature dependence of the elastic 
moduli cii and £33 at the phase transition in MnSi (the data 
are shown with offsets for better viewing). 

To avoid confusion we use the tilde sign to mark the elastic 
moduli to distinguish them from the constants of real 
tetragonal symmetry. 


II. EXPERIMENTAL 

In the course of the ultrasound studies of MnSi, quite 
a number of runs were performed using digital pulse-echo 


techniquesii. A disc shape single crystal sample of MnSi 
of '^15 mm in diameter and 2.15 mm in thickness with 
an orientation along the (100) plane was used. 

Two sample setups were employed, where the propa¬ 
gation vector k of the ultrasound wave was parallel or 
perpendicular to the direction of magnetic field. The 
propagation vector was always perpendicular to the (100) 
plane. The magnetic field was directed parallel or per¬ 
pendicular to the (100) plane. As a result, the demagne¬ 
tization factors for both setups appeared to be quite dif¬ 
ferent with the smaller value in the perpendicular setup, 
where the magnetic field was perpendicular to the prop¬ 
agation vector. 

Details of the sample preparation were described in^. 
The LiNbOs transducers were bonded to the sample with 
silicon grease. The temperature was measured using a 
calibrated Cernox sensor with an accuracy of 0.02 K. 

A sinusoidal pulse of ^40 MHz was sent to the trans¬ 
ducer that excited a sound wave inside the sample, which 
experienced multiple reflections. The sound travel time 
was determined by performing a cross-correlation be¬ 
tween two selected reflections. The speed of sound and 
elastic constants are then calculated using the known 
thickness and density of the samples and the relationship 
Cij = pV'^. The precision of the sound velocity determi¬ 
nations is no worse than one part in ^ 10® whereas the 
absolute accuracy is about 0.5%. The sound attenuation 
is estimated from an analysis of the amplitudes of two 
successive reflections as 20 /Lo^o 5 io(Ai/A 2 )- 


A. Longitudinal wave experiments 

The experimental data are presented in Figsl^HSl The 
temperature dependence of the elastic moduli cn, C33 
and corresponding data on the sound attenuation as a 
function of applied magnetic fields are shown in Figs. [2l 
[3] and m Note that application of a magnetic field to a 
cubic crystal makes it anisotropic. Accordingly we treat 
MnSi in a magnetic field directed along [001] as a tetrag¬ 
onal crystal (see also^^). However, to avoid confusion we 
prefer to use the tilde sign to mark the corresponding 
elastic moduli to distinguish them from the constants of 
real tetragonal symmetry. 

Fig. |S] displays the field dependence of cn and C33 at 
various temperatures. Distinct discontinuities of different 
signs at of the border of the conical and skyrmion phases 
most probably suggest a first order phase transitions^. 
It should also be pointed out that the elastic anisotropy, 
shown in Fig. [SI generally agrees with the data of RefSS, 
and corresponds to the uniaxial symmetry induced by 
magnetic field. 

The magnetic phase diagram of MnSi obtained as a 
result of these measurements is demonstrated in Fig. [51 
Note that the values of magnetic fields corresponding to 
some characteristic features of the diagram obtained in 
the two configurations, with magnetic field either parallel 
or perpendicular to the direction of sound wave propa- 
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FIG. 3: (Color online) Absorption of longitndinal sonnd waves 
at the phase transition in MnSi as fnnctions of magnetic field 
and temperature, k, u, H || [001] (the data are shown with 
offsets for better viewing). The numbers in the plot stand for 
values of magnetic field in Tesla. 

gation, are different. However, this situation is obviously 
connected to a difference in demagnetization factors aris¬ 
ing due to the disc shape of the sample, when the mag¬ 
netic field is directed parallel or perpendicular to the disk 
planei^. 

As is seen in Figs. [ 2 ] and 0 ] the magnetic phase transi¬ 
tion in MnSi in zero magnetic field is signified by a quasi 
discontinuity in cn, which suggests a first order character 
of the transformation. 

It is seen from Figs. [2] and 0 ] that the quasi discontinu¬ 
ities Acii and AC33 initially strongly decrease with mag¬ 
netic field, practically disappear in the range of ~0.1-0.2 
(~0.1-0.3) T, then recovers at ~0.2-0.45 (~0.3-0.5) T and 
again vanished at higher magnetic field. The numbers 
correspond to the measurements performed in the paral¬ 
lel and perpendicular setups, respectively. The method 
to evaluate the amplitude Acn and AC33 is illustrated in 
the inset of Fig. 0 ] 

This kind of evaluation of Acn and AC33 does not 
account for anomalous behavior different from quasi dis¬ 
continuities that also can be seen in Fig .01 Note that the 
region with almost zero values of Acn and AC33 in the 
interval ~0.12-0.26 (~0.12-0.2) T closely corresponds to 
the extent of the skyrmion phase along the magnetic to 
paramagnetic transition (Fig. 0 ]). 

This implies that the cn and C33 elastic constants 




Applied field (T) 


FIG. 4 : (Golor online) a) Amplitudes of cn and C33 disconti¬ 
nuities at the phase transition in MnSi as functions of mag¬ 
netic field. Inset in the upper right corner shows a way of 
determining the amplitudes, b) Partial amplitudes of absorp¬ 
tion of longitudinal sound waves at the phase transition in 
MnSi as functions of magnetic field. Inset in the upper left 
corner shows a way of determining the amplitudes. The per¬ 
pendicular setup k L H corresponds to the smaller demagne¬ 
tization factor. The difference between both sets of data a) 
and b) in Fig. 0 ] reflect mostly the difference in demagnetiza¬ 
tion factors for two positions of our disc like sample in regards 
to magnetic field. That was seen already in Fig. [2] and would 
also be seen in all subsequent figures. 



Applied field (T) 


FIG. 5: (Color online) Magnetic field dependence of the longi¬ 
tudinal elastic moduli at different temperatures. Temperature 
values for each isotherm starting from the top to the bottom 
are shown in the right lower corners of a) and b). 


are almost continuous through the transition from the 
skyrmion to the paramagnetic phase. Another notable 
feature of Fig. 0 ] are the maxima of Acii and the atten¬ 
uation at ~0.4(0.3) T, which correspond to the position 
of the announced tricritical point^^. The data on the 
sound attenuation (Figs. 01 Hb) confirm a separation of 
the boundary between paramagnetic and magnetically 
ordered phases into three distinct regions, which can be 
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FIG. 6 : (Color online) Magnetic phase diagram of MnSi ac¬ 
cording to the present studies. Closed circles-discontinuities 
of Cii and C 33 and maxima sound adsorption. Open circles- 
minima in cii and C33, Diamond-helical=conical transition, 
sf-strongly fluctuating region. 


named as helical, skyrmionic, and strongly fluctuating. It 
needs to be emphasized that the positions of the minima 
in cii(T) and C 33 (T) shift with magnetic field (Fig. [^l in 
such a way that they meet the phase transition line at 
-0.4 (0.3) T (Fig. HD. 


B. Shear wave experiments 



FIG. 7: (Color online) Shear modulus C44, measured in the 
( 001 ) plane with polarization u along the [ 100 ] and [ 010 ] di¬ 
rections. Here Aca = cu — Cii{H — 0, 27.45A'). The temper¬ 
atures of the isotherms, from top to bottom, correspond to T 
in K: 28.07, 27.85, 27.64, 27.45. 



Applied field (T) 


FIG. 8 : (Color online) Shear modulus C44, measured in the 
(001) plane with polarization u at an angle of 30 and 120 
degrees to the [100] direction. Here Aca = cu — ca {H = 
0,27.14A'). The temperatures of the isotherms, from top to 
bottom, correspond to T in K: 28.06, 27.84, 27.64, 27.44, 
27.24, 27.14. 


Experimental results on the propagation of the shear 
wave in MnSi are displayed in Figs. [THE Fig. [7] shows 
the shear elastic response of MnSi with both the propa¬ 
gation vector k and magnetic field H directed along [001], 
whereas polarization vector u is directed along [010] and 
[100]. 

One should remember that in the case of H = 0 this re¬ 
sponse corresponds to the C44 elastic constant even at an 
arbitrary direction of polarization u. In the current case 
of 7 ^ 0 a well-expressed anisotropy is quite evident. 

We practically cannot see the skyrmions at the indi¬ 
cated directions of polarization (Fig. [7D. Surprisingly, 
however, at some special positions skyrmion features are 
clearly exposed (Fig. ED. In contrast, in the perpendicu¬ 
lar setup of H II [001] and u [| [100] or [010] the skyrmions 
are clearly seen, see Fig. EJ 


III. DISCUSSION 

As was pointed out in the Introduction the heat ca¬ 
pacity, thermal expansion coefficient, and the tempera¬ 
ture coefficient of the electrical resistivity of MnSi at zero 
magnetic field are characterized by sharp peaks situated 
on the low temperature side of shallow maxima or min- 
imai^ii. The sharp peaks are considered now to be revela¬ 
tions of first order phase transitions, whereas the shallow 
maxima or minima arise as a result of intense helical fluc- 
tuations^id^. The peaks in the above-mentioned quanti¬ 
ties are the result of differentiations of step like functions, 
which describe the behavior of enthalpy, volume, and re¬ 
sistivity at first order phase transitions. At the same 
time the sound wave propagations directly probe values 
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FIG. 9: (Color online) Shear moduli £44 and cee measured 
in the (100) plane. Here Acu = cu — Cii{H = 0,27.9411'). 
The temperatures of the isotherms, from top to bottom, cor¬ 
respond to T in K: 28.85, 28.63, 28.43, 28.33, 28.24, 28.04, 
27.94. 


of elastic moduli and immediately give values of any dis¬ 
continuities in elastic properties. Therefore, the behavior 
of elastic moduli inferred from ultrasound studies better 
characterizes the nature of phase transitions than heat 
capacity or thermal expansions measurements. The lat¬ 
ter can be misleading when, for instance, a first order 
phase transition smears out by impurities or mechanical 
stresses. 

In the present case (see Figs. [5] and 0]) the elastic modu¬ 
lus cii is practically discontinuous with a positive change 
cii,p > ( p and h stand for paramagnetic and heli¬ 

cal) through the phase transition at zero magnetic field. 
The quasi discontinuities Acn decrease with magnetic 
field and become nearly zero across the phase transi¬ 
tion from the skyrmion to the paramagnetic phase, in 
the range ^0.12-0.26 (~0.12-0.2) T. 

A reasonable explanation of the observed facts would 
be that the paramagnetic phase coexisting with the 
skyrmion phase is dominated by skyrmion fluctuations. 
This implies that a skyrmion liquid with a longitudinal 
elasticity not much different from that of the skyrmion 
crystal may arise in MnSi on application of a magnetic 
field. Actually this could be verified in neutron scatter¬ 
ing experiments (strong chiral skyrmion like fluctuations 
were detected slightly above a.t H = 0 [9,19]). An¬ 
other explanation could be connected with an inefficiency 
of the Brazovski mechanism of appearance of a first or¬ 
der phase transition in case of the skyrmion-paramagnet 
transformationiS,. 

Another puzzle can be seen in Fig. 01 As one can see, 
the unexpected maximum of Acn and AC 33 , accompa¬ 
nied by an enhancement of the sound attenuation, is lo¬ 
cated in the vicinity of 0.4 (0.3) T, where the minimum 
in cii(T) and AC 33 approaches the phase transition line. 
This minimum originates from the strong helical fluctu¬ 
ations, which are believed to induce a first order phase 
transition in MnSi at least at zero magnetic field. Do 


these fluctuations help to recover first order features of 
the phase transition? Probably yes, because these fea¬ 
tures along with the sound attenuation decay rapidly in 
high magnetic fields, when the ’’helical” fluctuations are 
finally suppressed. (Note that the minima in cii(T) at 
B > 0.35 T are connected with a transition from the spin 
polarized paramagnetic state to a paramagnet.). In any 
case, a fluctuation mechanism of a first order phase tran¬ 
sition should stop working at low temperatures, when the 
thermal fluctuations die out. So if one believes that the 
phase transition in MnSi is first order at zero magnetic 
field then there should exist a corresponding tricritical 
point, where a first order transition meets a second order 
one. An existence of a tricritical point on the transi¬ 
tion line in MnSi at ~0.4 T and ^28 K was claimed in 
Repl®. Our ultrasound study cannot support this con¬ 
clusion. Indeed, an elastic modulus cn is expected to 
diverge at a tricritical point, which is not observed in the 
present investigation (see Fig. 

Taking a broader view of the situation one may come 
to the conclusion that the magnetic phase transition in 
MnSi is weak first order probably in the whole range of 
magnetic fields studied. Taking into account a first order 
nature of the phase line between conical and skyrmion 
phases, one may conclude that the crossing points 1,2 
(Fig. [T]:) are usual triple points, which are not accompa¬ 
nied by an anomalous behavior of elastic properties. The 
quasi discontinuities of cn and C33 vary from small to tiny 
under an influence of magnetic fields. Their recovery and 
the enhanced sound absorption occur at the intersection 
of the phase transition line and the line of minima in cn 
and C33. The latter is a region of strong spin fluctua¬ 
tions, reflected in maxima of heat capacity, etc. Despite 
the fact that this region does not carry features of a real 
phase transition, the existing powerful fluctuations at the 
minima of cn and C33 make the mentioned intersection 
point similar to a critical end point, where a second order 
phase transition meets a first order one. The best exam¬ 
ple of this situation is the crossing of the A-transition line 
in liquid helium and its melting curve, which is accom¬ 
panied by anomalies in the thermodynamic properties of 
the melting transition. 

One more enigma of the present study is the discovery 
of anisotropic behavior of the shear modulus (see Figs. [71 
B- The anisotropy was spotted at different positions of 
the polarization vector u in the ( 001 ) plane in the parallel 
setup. An evidence of the helical-conical transition at 
^0.1 T can be seen only at some selected positions of 
vector u. 

Another surprising detail is that the skyrmion phase 
does not reveal itself at the regular directions of the po¬ 
larization vector, along [100] or [010]. Instead, it can be 
clearly observed at angles of '^30 or 120 degrees for these 
directions, which certainly suggests a connection with the 
triangular packing. A skyrmion crystal is believed to be¬ 
long to this class of lattice. But this association can 
not be verified in the case of true hexagonal symmetry. 
However, the observed anisotropy is, as a matter of fact. 
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small. The difference of the shear moduli of the skyrmion 
phase measured in two directions of polarization vector is 
only ~ 2 X 10“^ % (see Fig. [ 8 ]). Hence, we may conclude 
that the hexagonal symmetry of the skyrmion crystal is 
slightly distorted by the magnetic field inhomogeneity or 
structural defects, which can easily be missed in neutron 
studies. On the other hand, the measurements of the 
shear moduli in the perpendicular setup do not show any 
noticeable anomalies. 

Getting to a discussion of the magnetic phase diagram 
of MnSi, we should note that the spotted scale difference 
in figures entitled as a and b is explained by a difference 
in the demagnetization factor. Indeed, in our sample the 
demagnetization factors differ by a factor of five for the 
parallel and perpendicular setups. At the same time, it 
needs to be emphasized that the skyrmion domain in case 
of the ’’perpendicular” setup with the smaller demagne¬ 
tization factor has a reduced temperature range, which 
suggests that the magnetic field inhomogeneity plays an 
important role in the skyrmion occurrence. 

IV. CONCLUSION 

The sound speed and attenuation were measured in a 
MnSi single crystal in the temperature range of 2 - 40 K 
and magnetic fields to 7 Tesla. 

The magnetic phase transition in MnSi in zero mag¬ 
netic field is signified by a quasi discontinuity in the cn 
elastic constant, which suggests a first order character of 
the transformation. The quasi discontinuities of cn and 
C 33 initially strongly decrease with magnetic field, prac¬ 
tically disappear in the range of ^ 0.1 — 02/0.3 T, then 
recover at ~ 0.2/0.3 — 0.45/0.5 T and again vanish at 
higher magnetic field. Note that the region with almost 
zero values of cn and C33 in the interval ~0.1 —0.2/0.3T 
closely corresponds to the extent of the skyrmion phase 
along the magnetic to paramagnetic transition. This im¬ 
plies that the cn and C33 elastic constants are almost con¬ 


tinuous through the transition from the skyrmion to the 
paramagnetic phase. This result probably questions the 
applicability of the Brazovsky theory for the skyrmion- 
paramagnet phase transitioiii^. The recovery of the dis¬ 
continuities of cii and C33 and the enhanced sound ab¬ 
sorption occur at the intersection of the phase transition 
line and the line of minima in cn and C33. The latter is a 
region of strong spin fluctuations, reflected in maxima of 
the heat capacity. Despite the fact that this region does 
not carry features of a real phase transition, the existing 
powerful fluctuations at the minima of Cn and C33 make 
the mentioned intersection point similar to a critical end 
point, where a second order phase transition meets a first 
order one. 

The spotted magnetic scale difference in the figures 
entitled a and b results from differences in the demag¬ 
netization factors. The skyrmion domain in case of the 
’’perpendicular” setup with the smaller demagnetization 
factor has a reduced temperature range, which suggests 
that the magnetic field inhomogeneity plays an impor¬ 
tant role in the skyrmion occurrence and, hence, opens a 
way of skyrmion manipulation (see also^^^. 

The elastic uniaxial symmetry of MnSi (Fig. [5]), in¬ 
duced by a magnetic field is broken by the magnetic field 
inhomogeneity or structural defects that also distort the 
hexagonal symmetry of skyrmion crystal, as evidenced by 
the anisotropy of the shear modulus in the (001) plane 
(Fig. [a ID). 
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